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Abstract. A multivariate polynomial is stable if it is nonvanishing whenever 
all variables have positive imaginary parts. We classify all linear partial differ- 
ential operators in the Weyl algebra An that preserve stability. An important 
tool that we develop in the process is the higher dimensional generalization 
of Polya-Schur's notion of multiplier sequence. We characterize all multivari- 
ate multiplier sequences as well as those of finite order. Next, we establish a 
multivariate extension of the Cauchy-Poincare interlacing theorem and prove 
a natural analog of the Lax conjecture for real stable polynomials in two vari- 
ables. Using the latter we describe all operators in that preserve univariate 
hyperbolic polynomials by means of determinants and homogenized symbols. 
Our methods also yield homotopical properties for symbols of linear stability 
preservers and a duality theorem showing that an operator in An preserves 
stability if and only if its Fischer-Fock adjoint does. These are powerful mul- 
tivariate extensions of the classical Hermite-Poulain-Jensen theorem, Polya's 
curve theorem and Schur-Malo-Szego composition theorems. Examples and 
applications to strict stability preservers are also discussed. 
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1. Introduction and main results 

In their seminal 1914 paper |49j Polya and Schur characterized all linear oper- 
ators that are diagonal in the standard monomial basis of C[z] and preserve the 
set of polynomials with only real zeros. Polynomials of this type and linear trans- 
formations preserving them are of central interest in e.g. entire function theory 
[17|, I37j : it is for instance well known that the Riemann Hypothesis is equivalent to 
saying that ^(^ + itj may be approximated by real zero polynomials uniformly on 
compact sets, where ^ denotes Riemann's xi-function. 

Polya-Schur's result generated a vast literature on this subject and related topics, 
see [1^ and references therein. Nevertheless, complete solutions to the fundamental 
problems of describing all linear operators preserving the set of real zero polynomials 
or, more generally, the set of polynomials with zero locus in a prescribed region 
51 C C, are yet to be found. Although many special cases and variations of these 
problems have been intensely studied for more than a century, to the best of our 
knowledge they have been stated in the above general (and explicit) form only 
recently by Craven-Csordas [17] and Csordas [T^, see also [Tj Ffl fTOl [T5] . 

This paper is part of a series [3 HI [ill El IE] devoted to these questions, 
their natural multivariate extensions and applications to geometric function theory, 
matrix theory, probability theory, combinatorics, and statistical mechanics. Here 
we 

(1) classify all linear partial differential operators in the n-th Weyl algebra that 
preserve stable, respectively real stable polynomials in n variables; 

(2) obtain a Lax type determinantal representation for linear operators in the 
first Weyl algebra (n = 1) preserving real zero polynomials and a charac- 
terization in terms of their homogenized symbols; 

(3) prove higher dimensional versions of Polya-Schur's theorem; 

(4) apply (l)-(3) to establish a Fischer-Fock duality for stability preservers, ho- 
motopic properties of their symbols and geometric interpretations extending 
Polya's "curve theorem" for all n > 1, stable multivariate generalizations of 
the Cauchy-Poincare interlacing theorem, Schur-Malo-Szego type convolu- 
tion theorems in higher dimensions, and necessary and sufficient conditions 
for strict (real) stability preserving in one or several variables. 

A nonzero univariate real polynomial with only real zeros is called hyperbolic 
while / e C[z] is called stable if f{z) ^ for all z G C with 3m{z) > 0. Hence a 
univariate real polynomial is stable if and only if it is hyperbolic. These classical 
concepts have several natural extensions to multivariate polynomials, see, e.g., four 
different definitions in [33j- Below we concentrate on the most general notion: 

Definition 1.1. A polynomial / S C[zi, . . . , Zn] is stable if /(zi, . . . , z„) for 

all 71-tuples (zi, . . . , Zn) G C" with 3m{zj) > 0, 1 < j < n. If in addition / has real 
coefficients it will be referred to as real stable. 

Clearly, / is stable (respectively, real stable) if and only if for all a e M" and v g 
R" the univariate polynomial /(a -I- vt) G C[t] is stable (respectively, hyperbolic), 
see Lemma \n\ in Section [51 In what follows we denote by Hn{C), respectively 
Hn(M), the set of stable, respectively real stable polynomials in n variables. 

Another fundamental extension of the notion of real-rootedness to higher dimen- 
sions stems from PDF theory. Namely, a homogeneous polynomial p € K[zi, . . . , Zn] 
is said to be (Carding) hyperbolic with respect to a given vector u S R" ]ip{v) ^ 
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and for all vectors a G M" the univariate polynomial p{a + vt) e M[t] has only real 
zeros. For background on (multivariate homogeneous) hyperbolic polynomials one 
may consult, e.g., In Section [^?^ we prove the fohowing result describing 

the relation between real stable and hyperbolic polynomials. 

Proposition 1.1. Let f E M.[zi, . . . , z„] be of degree d and let Jh G K[zi, . . . , Zn+i] 
he the (unique) homogeneous polynomial of degree d such that ///(zi, . . . , z„, 1) = 
/(zi, . . . , Zn). Then f € 7i„(R) if and only if fn is hyperbolic with respect to every 
vector V G R"+^ such that Vn+i = and Vi > 0, 1 < i < n. 

It is worth mentioning that real stable multivariate polynomials appear already 
in Theorem 1 of the foundational article [26] by Garding and that stable multivariate 
entire functions can be found in Chap. IX of Levin's book [37]. 

Let An [C] be the Weyl algebra of all finite order linear differential operators with 
polynomial coefficients on C[zi, . . . , z„]. Recall the standard multi- index notation 

= z^' = di'---dn", where z = (zi,...,z„), a = (ai,...,a„) G N" 

and di = d/dzi for 1 < i < n. Then each operator T S ytn[C] may be (uniquely) 
represented as 



where Oafs € C is nonzero only for a finite number of pairs (a, (3). Let further .4„[R] 
be the set of all T G .4„[C] with Oafj G K for all a, /3 G N". A nonzero differential 
operator T G An[C\ is called stability preserving if T : H„(C) H„(C) U {0} and 
it is said to be real stability preserving if T : 7Y„(M) H„(M) U {0}. 

Given T of the form ([!]) define its symbol Ft{z,w) to be the polynomial in 
C[zi, . . . ,z„,-u;i, . . . ,Wn] given by Ft{z,w) = J2a,t3 "'a(3Z°'w'^ . 

The first main results of this paper are the following characterizations of the 
multiplicative submonoids .4„(C) C .A„[C] and -4,i(M) C .4n [M] consisting of all 
stability preservers and real stability preservers, respectively. 



Theorem 1.2. Let T G Ai[C]. Then T G Ai(C) if and only if Ft{z,-w) G 

H2„(C). 

Theorem 1.3. Let T G Ai[R]- Then T G An(M) if and only if Ft{z,-w) G 

W2„(M). 



It is interesting to note that Theorems 11.21 and 11.31 essentially assert that finite 
order stability (respectively, real stability) preservers in n variables are generated 
by stable (respectively, real stable) polynomials in 2n variables via the symbol map. 
Geometric interpretations of these statements in terms of symbol surfaces are given 
in Section 

To prove the above theorems we need to generalize a large number of notions 
and results for univariate stable and hyperbolic polynomials to the multivariate 
case. Let ai < a2 < • ■ ■ < «„ and A < /32 < • ■ ■ ^ /^m be the zeros (counted 
with multiplicities) of two given polynomials f,gG Hi (M) with deg f — n and 
degg = m. We say that these zeros interlace if they can be ordered so that either 
0(1 < Pi < ct2 < /32 < • ■ ■ or pi < ai < P2 < 012 < • • • , in which case one 
clearly must have \n — m\ < 1. Note that by our convention, the zeros of any two 
polynomials of degree or 1 interlace. It is not difficult to show that if the zeros of 
/ and g interlace then the Wronskian W[f, g] := f'g — fg' is either nonnegative or 
nonpositive on the whole real axis M, see, e.g., [51]. In the case when g] < 




(1) 



a,/3eN' 
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we say that / and g are in proper position, denoted f <^ g. For technical reasons 
we also say that the zeros of the polynomial interlace the zeros of any (nonzero) 
hyperbolic polynomial and write <C / and f <ti 0. Note that if /, g are (nonzero) 
hyperbolic polynomials such that f <^ g and g <ti f then / and g must be constant 
multiples of each other, that is, = 0. 

The following theorem is a version of the classical Hermite-Biehler theorem [5T] . 

Theorem 1.4 ((Hermite-Biehler theorem)). Let h := f + ig (z C[z], where f,gG 
R[z]. Then h € 7ii(C) if and only if g <^ f . 

The Hermite-Biehler theorem gives an indication about how one should gener- 
alize the concept of interlacing to higher dimensions: 

Definition 1.2. Two polynomials /, <? G ^[zi, ■ ■ ■ ,Zn] are in proper position, de- 
noted / < 5, if .9 + e Wn(C). 

Equivalently, / and g are in proper position if and only if for all a G M" and 
V € K" the univariate polynomials f(a + vt), g{a + vt) G are in proper position. 
It also follows that f,gG 7i„(R) U {0} whenever f g, see Corollarv l2.4l in Section 



The next (also classical) result is often attributed to Obreschkoff [45] and some- 
times referred to as the Hermite-Kakeya-Obreschkoff theorem [51] . 

Theorem 1.5 ((Obreschkoff theorem)). Let /, g e R[z]. Then af + f3g e HiiR) U 
{0} for all a, /? G K i/ and only if either f 'ti g, g ^ f or f = g = 0. 

We extend this theorem to polynomials in several variables as follows. 

Theorem 1.6. Let f,g £ R[zi, . . . , z„]. Then af + Pg e Hn{^) U {0} for all 
a, [3 gR if and only if either f g , g f or f ^ g = 0. 

Recall that an infinite sequence of real numbers A : N — > M is called a mul- 
tiplier sequence (of the first kind) if the associated linear operator T on C[z] 
defined by T{z") — X{n)z", for all n e N, is a hyperbolicity preserver, i.e., 
T : Hi(M) Hi(M)U{0}. Any linear operator T : C[zi, . . . , z„] ^ C[zi, . . . , z„] can 
be represented as a formal power series in d with polynomial coefficients. Indeed, 
this may be proved either by induction or by invoking Peetre's abstract characteri- 
zation of differential operators 37] . Note also that in general a multiplier sequence 
is represented by an infinite order differential operator with polynomial coefficients. 

In [49| Polya and Schur gave the following characterization of multiplier se- 
quences of the first kind. 

Theorem 1.7 ((Polya-Schur theorem)). Let X : N ^ R be a sequence of real 
numbers and T : R[z] —^ R[z] be the corresponding (diagonal) linear operator. 
Define ^{z) to be the formal power series 



The following assertions are equivalent: 

(i) X is a multiplier sequence, 

(ii) ^{z) defines an entire function which is the limit, uniformly on compact 
sets, of polynomials with only real zeros of the same sign, 



m 




k=0 
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(iii) Either $(z) or $(— z) is an entire function that can be written as 

oo 

where n € N, C £ M, a, > for all k Cz N and X]fc°=i < oo, 

(iv) For all nonnegative integers n the polynomial T[{1 + z)"] is hyperbolic with 
all zeros of the same sign. 

We introduce a natural higher dimensional analog of the notion of multiplier 
sequence and completely characterize all multivariate multiplier sequences as well 
as those that can be represented as finite order differential operators. For this we 
need the following notation. Given an integer n > 1 and a, P € N" we write a < (3 
if Ui < Pi iov 1 < i < n. Let further \a\ = |ai| + • • • + |a„|, = a^^ ■ ■ ■ af^" , 
a! = ai! • • • q;„!, = if a ^ /3 and = /?!/(/? — a)! otherwise. 

Definition 1.3. A function A : N" ^ M is a (multivariate) multiplier sequence if 
the corresponding (diagonal) linear operator T defined by T(z") = X{a)z°', for all 
a € N", is a real stability preserver, that is, T : 7^„(M) -> 7i„(M) U {0}. 

The following theorem completely describes multivariate multiplier sequences. 

Theorem 1.8. Consider an arbitrary map X : N" —>■ M. Then X is a multivariate 
multiplier sequence if and only if there exist usual (univariate) multiplier sequences 
Xi : N ^ M, for all 1 < i < n, such that 

X{a) = Ai(q!i)A2(q!2) ■ • • Xn{an), for all a = (ai, . . . , «„) e N", 

and either A(a)A(/3) > for all a,/3 G N", or (-l)l"l+l^l A(a)A(/3) > for all 
a,pe W. 

Note that Theorem 11.81 is a negative result, since it shows that the only mul- 
tivariate multiplier sequences are those that one would expect exist: products of 
univariate ones. 

We next characterize all multiplier sequences that are finite order differential 
operators, i.e., those whose symbols are (finite degree) polynomials: 

Theorem 1.9. Given a map X : N" R, let T be the corresponding (diagonal) 
linear operator. Then T € .A„(R) if and only if T has a symbol Ft{z,w) of the 
form 

Ft{z, w) = fl{ziWl)f2{z2W2) ■ ■ ■ fn{ZnWn), 

where fi{t) is a polynomial with only real and nonpositive zeros for all 1 < i < n. 

Remark 1.1. Note that Theorem 11.91 combined with well-known properties of uni- 
variate multiplier sequences (cf. Lemma 13.11 below) implies in particular that if 
A : N" — > M is a finite order multivariate multiplier sequence, then there exists 
7 G N" such that A(a) = for a < 7 and either A(a) > for all a > 7 or A(a) < 
for all a > 7. Note also that for n = 1 Theorem ll.QI gives an alternative description 
of finite order multiplier sequences that complements Polya-Schur's Theorem 11.71 

Our next result is a vast generalization of the following classical theorem |5T] . 

Theorem 1.10 ((Hermite-Poulain- Jensen theorem)). Let p{z) — X]fc=o '^^-^'^ ^ 
K[z] be nonzero and let T = X]fe=o CLkd*' / dz^ € ^i[R]. Then T € yli(R) if and only 

ifpeHiiR). 
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The natural setting for our extension is the Fischer-Fock space !Fn [211 1221 1231 
I24j , also called the Bargmann- Segal space [3l |4l [55] or the Newman- Shapiro space 
[321 1331 [331 [5^ , which is the Hilbert space of holomorphic functions / on C" such 
that 



Here a{a)z°' is the Taylor expansion of /. The inner product in Tn is given by 



where 5 is the Kronecker delta and denotes the z-th standard generator of the 



lattice Z". Hence, if T = Y^^.p a^pz^'d'^ € AJC] then 

{T{f),g) = ^a„^(z"9/'/,5) = <^Mf , d'^ 9) = (/, J] o^^^S'^s). 



Therefore, the formal Fischer- Fock dual (or adjoint) operator of T is given by T* = 
/3^/3^z"9'', cf. ;43J. Note that for 1 < j < n the dual of di is the operator given 
by multiplication with Zi and that diagonal operators (in the standard monomial 
basis) are self-dual. In particular, if T is a multiplier sequence then T* = T. 

Remark 1.2. The Fischer-Fock space J-n was used by Dirac to define second quan- 
tization [20] and its inner product has since been rediscovered in various contexts, 
e.g. in number theory where the corresponding norm is known as the Bombieri 
norm [5ll53j. Further important properties of J-n such as its (Bergman- Aronszajn) 
reproducing kernel and the Newman-Shapiro Isometry Theorem may be found in 
[42 l l43 l l44] . We should also point out that in e.g. P- module theory and microlocal 
Fourier analysis [40j one usually works with the inner product on J^n defined by 
{f{z), g{z))d — {f{iz),g{iz)), where (•,•) is as in ([2]). Note that the dual operator 
of di with respect to (•, ■)d is the operator given by multiplication with — z^. 

In Section [3131 we give a geometric interpretation and proof of the fact that the 
duality map with respect to the above scalar product preserves both ^„(C) and 
An (K) ■ More precisely, from Theorems I1.2H1.3I we deduce the following natural 
property: 

Theorem 1.11 ((Duality theorem)). Let T G An[C]. Then T £ ^„(C) if and only 
ifT* e An{C). Similarly, ifT € AnW\ then T e A»(K) if and only ifT* G An(M)- 

We conclude this introduction with a series of examples of real stable polynomials 
and various applications of our results. Further interesting examples of multi-affine 
stable and real stable polynomials can be found in e.g. [9l [TT l [M l I16j. 

Proposition 1.12. Let Ai, . . . , An be positive semidefinite mxm matrices and let 
B be a complex Hermitian mxm matrix. Then the polynomial 






(3) 



is either real stable or identically zero. 
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A proof of the above proposition is given in Section [6l Using this result and 
Theorem 11.21 we obtain a multidimensional generalization of the Cauchy-Poincare 
interlacing theorem: see Theorem 16.21 in Section l6.ll 

The Lax conjecture j35| for (Garding) hyperbolic polynomials in three variables 
has recently been settled by Lewis, Parillo and Ramana [38]. Their proof relies on 
the results of Helton and Vinnikov 30J. Applications of these results to e.g. hy- 
perbolic programming and convex optimization may be found in [52j . In Section 
16.21 we prove the following converse to Proposition 11.121 in the case n = 2 and thus 
establish a natural analog of the Lax conjecture for real stable polynomials in two 
variables. 

Theorem 1.13. Any real stable polynomial in two variables x,y can be written as 
±det{xA + yB + C) where A and B are positive semidefinite matrices and C is a 
symmetric matrix of the same order. 

Remark 1.3. A characterization of real stable polynomials in an arbitrary number 
of variables has recently been obtained in [14 . 

Combining Theorem 1 1 . 1 31 with Theorem 1 1 . 31 we get two new descriptions of finite 
order linear preservers of hyperbolicity (i.e., univariate real stability), namely a 
determinantal characterization and one in terms of homogenized operator symbols: 
see Theorems 16.81 and 16 .91 in Section [6731 

Further applications of our results include multivariate Schur-Malo-Szego com- 
position formulas and closure properties under the Weyl product of (real) stable 
polynomials fSection 14. 3p . a unified treatment of Polya type "curve theorems" as 
well as multivariate extensions fSection 14. 4p . and necessary and sufficient criteria 
for strict stability and strict real stability preservers (Section [5]). 

Brief excursion around the literature. The study of univariate stable polynomi- 
als was initiated by Hermite in the 1860's and continued by Laguerre, Maxwell, 
Routh, Hurwitz and many others in the second half of the XlX-th century. The 
contributions of the classical period are well summarized in [25l [501 [S] • Impor- 
tant results on stability of entire functions were obtained in the mid XX-th century 
by e.g. Krein, Pontryagin, Chebotarev, Levin [37]. Modern achievements in this 
area can be found in [46] and references therein. Much less seems to be known 
concerning multidimensional stability. In control theory one can name a series of 
papers by Kharitonov et al. [33j with numerous references to the earlier litera- 
ture on this topic. Another origin of interest to multivariate stable polynomials 
comes from an unexpected direction, namely the Lee- Yang theorem on ferromag- 
netic Ising models, the Heilmann-Lieb theorem for monomer-dimer systems and 
their various generalizations [5S1 [3S1 [35] ■ Combinatorial theory provides yet an- 
other rich source of stable polynomials as multivariate spanning tree polynomials 
and generating polynomials for various classes of matroids turn out to be stable 
(cf., e.g., [m [16]). Multivariate stable polynomials were recently used in [28] to 
generalize and reprove in a unified manner a number of classical conjectures, in- 
cluding the van der Waerden and Schrijver- Valiant conjectures, and in [S] to solve 
some long-standing conjectures of Johnson and Bapat in matrix theory. Further 
recent contributions include [8] , where a complete classification of linear preservers 
of univariate polynomials with only real zeros - and, more generally, of univariate 
polynomials with zeros only in a closed circular domain or on the boundary of such 
a domain - has been obtained, thus solving an old open problem going back to 
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Laguerre [34j and Polya-Schur [4^. Let us finally mention that real stable polyno- 
mials have also found remarkable applications in probability theory and interacting 
particle systems. Indeed, these polynomials were recently used in '11] to develop a 
theory of negative dependence for the class of strongly Rayleigh probability mea- 
sures, which contains several important examples such as uniform random spanning 
tree measures, fermionic/determinantal measures, balls- and-bins measures and dis- 
tributions for symmetric exclusion processes. 

Acknowledgements. It is a pleasure to thank the American Institute of Math- 
ematics for hosting the "Polya-Schur-Lax Workshop" on these themes in Spring 
2007. We would also like to thank an anonymous referee for several suggestions on 
improving the exposition of this paper. 

2. Basic properties and generalized Hermite-Kakeya-Obreschkoff 



The following criterion for (real) stability is an easy consequence of the defini- 
tions. 

Lemma 2.1. Let f e C[zi, . . . . z„]. Then f G 7Y„(C) if and only if f{a + vt) G 
Hi(C) for all a G K" and veRl- 

The next lemma extends the Hermite-Biehler theorem to the multivariate case 
and provides a useful alternative description of the proper position/ "interlacing" 
property for multivariate polynomials. 

Lemma 2.2. Let f,g G M[zi, . . . , z„] and let z„+i be a new indeterminate. Then 
f g if and only if g + Zn+i/ G "Hn+ii^)- Moreover, if f £ 7i„(M) then f ^ g if 
and only if 



whenever Jm{zi) > for all 1 < i < n. 

Proof. The "if" direction is obvious. Suppose that f g and that z„+i = a + ib, 
where a G R and b G K+. Then by Lemma l^TTl we have that f{a+vt) ^ g{a+vt) for 
all a G M" and v G M" . By Obreschkoff's theorem the zeros of g{a + vt) + af{a + vt) 
and bf{a + vt) interlace (both cannot be identically zero). Moreover, 

W{bfia + vt),g{a + vt) + af{a + vt)) = bW{f{a + vt)), g{a + vt)). 

Thus bf{a + vt) < g{a + vt) + af{a + vt) for aU a G M" and w G M!f:, which by 
Lemma [2. II implies that g + {a + ib)f G 7i„(C). But g + Zn+if clearly has real 
coefficients so 17 -I- Zn+if G 7i„+i(M). The final statement of the lemma is a simple 
consequence of the above arguments. □ 

Lemma 2.3. Suppose that fj G C[zi, . . . , z„] for all j G N are nonvanishing in an 
open set U C C" and that f is the limit, uniformly on compact sets, of the sequence 
{fj}jeN- Then f is either nonvanishing in U or it is identically equal to 0. 

Proof. The lemma follows from the multivariate version of Hurwitz' theorem on 
the continuity of zeros of analytic functions, see, e.g., |16j . □ 



Let /(zi, . . . ,z„) G Hn{C), a G M and A > 0. Then f{a + \zi, . . . ,z„) G H„(C). 
By letting A we have by Lemma [^31 that f{a, Z2, . . . , z„) G H,i_i(C) U {0}. 



Theorem 
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Corollary 2.4. For each n e N one has 

HniC) = {g + if:f,ge 7^„(R) U {0}, / « g}. 

Proof. The only novel part is that f,gE ?i„(M) U {0} whenever f g. This follows 
from Lemma [2.2l and Lemma [2.3l when we let z„+i tend to and cxd, respectively. □ 

We are ready to prove our multivariate Obreschkoff theorem. 

of Theorem[rB Suppose that f g. By Corollary [m we have g S '^„(M) U {0} 
so we can normalize and set /? = 1. By Lemma we have g + Zn+if G 7i„+i(R) C 
7i„+i(C), so by letting Zn+i = i + a with a e R we have g + af + if E 7i„(C), i.e., 
f -€1 g + af. From CoroUarv 12 . 41 again it follows that g + af E 7i„(R) U {0}, as was 
to be shown. 

To prove the converse statement suppose that we do not have f ~ g = 0. If 
af + Pg € n„{R) U {0} for aU a,P € R then (by Lemma [Q and Obreschkoff 's 
theorem) for all 7 € R" and v G R" we have either f{'^ + vt) <C 3(7 + vt) or 
f{j+vt) » g{j+vt). If both instances occur for different vectors, i.e., /(7i+wii) <C 
gili+vit) and 7(72 +f2i) > 5(72 +^^2*) for some 71,72 G R" and vi,V2 S R+, then 
by continuity arguments there exists r G [0, 1] such that /(7r + Vrt) <C gi'jr + Vrt) 
and /(7r + Wr^) ^ gilr + Vrt), whcrc 7t- r7i + (1 — t)72 G M" and u,- := 
TUi + (1 — t)v2 e M" . This means that /(7t + Vrt) and .g(7r + are constant 
multiples of each other, say /(7r + Vrt) — Xg{'-fr + v^t) for some A G K. By 
hypothesis we have h :— f — Xg E 7Y„(R) U {0} and h{'-fr + Vrt) = 0, in particular 
h{'yT + ivr) — 0. Since G M" it follows that /i = and / = \g. Consequently, if 
both instances occur we have f g for trivial reasons. Thus we may assume that 
only one of them occurs. But then the conclusion follows from Lemma [2 .11 □ 

Define 

HniCy = {/ G C[zi, . . . , z„] : /(zi, , . . . , z„) ^ if 3m{z,) < for 1 < ^ < n}. 

Clearly, f{z) E HniC) if and only if f{-z) E HniC)-. Hence by Corollary [23 and 
Lemma |2. II we have h := g + if E 7i,i(C)^ with f,gE R[zi, . . . , z„] if and only if 
5< /• 

Proposition 2.5. For any n E N the following holds: 

HniC) n HniC)- = CHniR) {cf : c E <CJ E -Hnim- 

Proof. Suppose that h = g+if E 7i„(C)n7Y„(C)" . By CoroUarvl^we have / < g 
and g ^ f . Hence for all a E R" and v E W\_ we also have f{a + vt) g{a + vt) 
and g{a + vt) <C f{a + vt). This means that /(a + vt) and g{a + f<) are constant 
multiples of each other, say f{a + vt) = \g{a + vt). By the multivariate Obreschkoff 
theorem we have that f — Xg E 7Y„(R) U {0}. Since (/ — Xg){a + vi) = we must 
have f -Xg = 0, i.e., /i = (1 + iX)g E C7i;„(M). □ 

Lemma 2.6. Let f E 7^„(R). Then the sets 

{g E Hn{R) ■ f « .9} and {g E H„(R) : / » ff} 

are nonnegative cones, i.e., they are closed under nonnegative linear combinations. 

Proof. Let / E 7i„(M) and suppose that f g and f <^ h. Then by Lemma [2T2l 
we have 3m{g{z) / f{z)) > and 3m(/i(z)//(z)) > whenever 3m{z) > 0. Hence if 
X,fJ.>0 and 3m{z) > then 3m((A(7(2;) + fih{z)) / f (zj) > and Lemma [221 yields 
f ^ Xg + fih. The other assertion follows similarly. □ 
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3. Classifications of multivariate multiplier sequences and finite 

ORDER ones 

3.1. Univariate and multivariate multiplier sequences. Let us first recall a 
few well-known properties of (usual) univariate multiplier sequences, see, e.g., |17) . 

Lemma 3.1. Let A : N — > M 6e a multiplier sequence. IfO<i<j<k are such 
that X{i)X{k) ^ then A(j) 7^ 0. Furthermore, either 
(i) all nonzero A(«) have the same sign, or 

ill) all nonzero entries of the sequence {(— l)'A(z)}i>o have the same sign. 
In what follows we denote the standard basis in M" by {ck : \ <k <n}. 

Remark 3.1. Suppose that a £ N", 1 < k < n, f{zk) := '^iLo'^i^l ^ and 
assume that A is a multivariate multiplier sequence. Then 

AT 

T{z^f{zk)) = ^ A(« + iek)a,zi e W„(K). 

Hence the function i 1— > \{a + ick) is a univariate multiplier sequence. 

The proofs of our characterizations of multivariate multiplier sequences and those 
of finite order build on a series of statements that we proceed to describe. 

Lemma 3.2. Let /(zi, 22) — Qoo + 0-01^2 + Oio^i + 011^122 G Z2] \ {0}- Then 
f E 7^2(1^) i/ imd only i/det(ay) < 0. 

Proof. Let a e M and denote hy A = (oij) the matrix of coefficients of /(zi,Z2). 
Clearly, /(zi, Z2) € ^2(1^) if and only if f{zi + a, Z2) G 7^2(1^) and /(zi, Z2 + a) G 
7^2(1^). We get the matrix corresponding to /(zi + «, Z2) by adding a times the last 
row oi A to the first row oi A, and we get the matrix corresponding to /(zi, Z2+a) by 
adding a times the last column of A to the first column of A. Since the determinant 
is preserved under such row and column operations we can assume that A has one 
of the following forms: 



/ floo 






aoiA 


/ floo 


aoi \ 


I 


ail ' 


Vaio 


i' 


I 


J 



Obviously, these matrices correspond to a polynomial /(zi, Z2) G 7^2(1^) if and only 
ifdet(ay)<0. □ 

Lemma 3.3. Let A : N" ^ M, 71 > 2, be a multivariate multiplier sequence and let 
7 G N" and I < i < j < n. Then 

A(7)A(7 + e, + gj) = A(7 + e,)A(7 + e^). 

Proof. Without loss of generality we may assume that i = 1 and j — 2. Let 
/(z) = z'^giz) G 7Y„(M), where .g(zi, Z2) = aoo + aoiZ2 + aiozi + aiiZiZ2 G 7Y2(K). 
It follows that the polynomial 

A(7)aoo + A(7 + e2)aoiZ2 + A(7 + ei)aioZi + A(7 + ei + e2)aiiZiZ2 

is stable or identically zero. By choosing A as 

(1 1) ^^^^ (1 -0 ' 

respectively, we get by Lemma [??^ that A(7)A(7 + ei +62) < A(7-|-ei)A(7-|-e2) and 
A(7)A(7 + ei + 62) > A(7 + ei)A(7 + 62), respectively, which proves the lemma. □ 
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Given a, /3 e N" with a < /3 set [a, /3] {7 e N" : a < 7 < /?}. 

Lemma 3.4. Let A : N" — > R fee a multivariate multiplier sequence. If a, (3 G N" 
are suc/i i/iaf A(a)A(/3) 7^ and 7 £ [a, f3] then A(7) 7^ 0. 

Proof. We use induction on ^ = |/3| — |a|, the length of the interval [a,/3]. The 
cases £ = and £ = 1 are clear. By Remark l3.1l and Lemma [3.11 the result is true 
in the univariate case. So we may assume that a and f3 differ in more than one 
coordinate, i.e., that a + ei, a + 62 G [a, /?]. 

If there exists an atom a + e^ € [a, (3] such that A(q; + Ci) =/= then by induction 
we have that A is nonzero in [a + e^, /?] . If a + Cj is another atom then a + e^ + Cj G 
[a + Ci, (3] so X{a + e^ + ej) ^ 0. Lemma [3.31 then gives that A(q; + ej) ^ 0. Thus, 
by induction, A is nonzero in [a + ej, P] for all a + Cj G [a, /?] and we are done. 

In order to get a contradiction we may assume by the above that A(q; + Ci) = 
for all a + Ci G [a, f3]. Let 7 G (a, (3] be a minimal element such that A(7) 7^ 0. If 
T is the (diagonal) linear operator associated to A then 

T(z"(l + z)''-") = A(a)z" + X{-f)z-' G 7i„(M). 

By Lemma [3?3l we have that X{a + ei + ej) = for all atoms a + Ci, a + ej G [a, /?]. 
By Remark 13.11 and Lemma 13.11 we also have A(a + me^) = for all to > 1 and 
atoms a + e; G [a,/?]. It follows that I7I — |q;| > 3. Now, if we set = t for all i 
then by the above we obtain that the polynomial 

A(a)tl"l +A(7)t'^' 

is hyperbolic in t, which is a contradiction since I7I — |a| > 3 and A(q;)A(7) 7^ 0. □ 

Lemma 3.5. Let A 

that X{a) 7^ 0. Then 



for all P > a. 

Proof We claim that A(7)A(7 + aei + feej) = A(7 + aei)A(7 + feej) for all 7 G N" and 
positive integers a and b. This follows easily by induction on a + fe using Lemma 13.31 
and Lemma 13.41 We proceed to prove the ^ by induction on the number r of 
nonzero entries oi 6 = (3 — a. The basis of induction, r < 1, is trivial. For r > 2 let 
j be an index such that 5j > and consider 7 — a + SjCj. By Lemma 15^ A (t ) 7^ 0. 
By induction 

XiPl^XiPl A(7) A(a+(/3,-7»)e.) A(a + (7, - a,)e,) 

A(a) A(7) ■ A(a) A(7) ' X{a) 

For i ^ j the corresponding factor is 

X{a + SiCi + SjCj) X{a) X{a + SiCi) 
X{a + 6jej) X{a) S{a) 
by the claim above. For i ~ j the corresponding factor is 

X{a + SjCj) X{a + djej) X^a + SiCj) 

X{a + SjCj) X{a) S{a) 
and the lemma follows by induction. 

□ 



N" ^ R fee a multivariate multiplier sequence and suppose 
X{a) li A(a) ^ ' 
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Figure 1 . Illustration of the induction step in Lemma 13.61 

If a = (ai, . . . , an) € N" and P = . . . , /3„) G N" we define two new vectors 
a V /3, a A /3 € N" by setting a V /? = (max(Q;i, . . . , max(a„, /?„)) and a A P = 
(min(ai, . . . , min(a„, /?„)). 

Lemma 3.6. Let A : N" W be a multivariate multiplier sequence and suppose 
that Xia)\{P) ^ 0. Then X{a y (}) ^ if and only if \{a A /?) ^ 0. 

Proof. If A(a)A(/3) ^ and A(q; A /3) 7^ then Lemma [3.51 and Lemma [3.11 imply 
that A(a V /?) ^ 0. Suppose that A(a)A(/3) ^ and A(a V /3) 7^ 0. We prove that 
A(q! a /?) 7^ by induction on |a — /?|. If ja — /3| = there is nothing to prove. Also, 
if a and (3 are comparable there is nothing to prove, so we may assume that there 
are indices i and j such that ai < (3i and [3j < aj . Since a<a + ei<a\/f3 and 
(3 < P + ej < a V /3 we have by Lemma \3A\ that X{a + ei)A(/3 + Cj) ^ 0. Consider 
the pairs (a + e^,/? + ej), {a + e^,/?) and (a,/3 + Cj). The distance between each 
of them is smaller than ja — /3|, they all have to join a V /3, and the meets are 
a A P + ei + Cj, a A P + Cj and a A P + Ci respectively, see Fig. [TJ By induction we 
have that X{a AP + ei)X(a AP + ej)X(a A /3 + + Cj) 7^ 0. By Lemma [3731 this gives 

X{aAp + e,)X{aAP + ej) 

X{a A P) — r-7 -; r 7^ (J, 

^ ' X{aAP + er + ej) ^ 

which is the desired conclusion. □ 

Recall that the support of a map A : N" ^ K is the set {a e N" : X{a) 7^ 0}. 

Lemma 3.7. Let X : N" — > M 6e o multivariate multiplier sequence. Then there 
exist univariate multiplier sequences Ai : N ^ R, 1 < * < J^, such that 

X{a) = Ai(q!i)A2(q!2) ■ • ■ A„(a„), a = (ai, . . . ,a„) G N". 

Proof. Let S be the support of A. By Lemma 13.51 and Remark 13.11 it suffices to 
prove that S has a unique minimal element. So far, by Lemma l3.4l and Lemma l3.6l 
we know that 5' is a disjoint union S = U™ of boxes Bi = !{ x ■ ■ ■ x where 

is an interval (possibly infinite) of nonnegative integers. Also, points in different 
boxes are incomparable. 

Suppose that S does not have a unique minimal element. We claim that there 
exists an interval [a, P] such that [a, P]r\S = {S, 7}, where S and 7 are in different 
boxes. We postpone the proof of this statement for a while and show first how it 
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leads to a contradiction. Let T be the (diagonal) linear operator associated to A. 
We then have 

r(z"(l + z)^-") = \{5)z^ + A(7)zT. 
Now |(5 — 7I > 3 since otherwise 5 and 7 would be comparable or we would have 
7 = (5 A 7 + ei and (5 = (5 A 7 + for some i and j . This is impossible by Lemma 
13.31 since 7 and 6 would then be in the same box. By assumption we have that 

A((5)z*-*^^ + A(7)z'^"^^^ e 7i„(R) 

so by setting all the variables in z^'^'^'^ equal to t and setting all the variables in 
^7-i5a7 gqyg^i —^-1 (which we may since z^~^'^'^ and 2;T-*^t contain no common 
variables) we obtain that 

A((5)tl*-^l ±A(7) eHi(M). 

This is a contradiction since — 7I > 3 and (5, 7 g S so A (5) A (7) ^ 0. 

It remains to prove the claim. Let d be the minimal distance between different 
boxes and suppose that i5 and 7 are two points that realize the minimal distance. 
If K € [(5 A 7, (5 V 7] then |k — (5| < d with equality only if k = 7 and |k — 7I < d with 
equality only li n — 5. It follows that [(5 A 7, (S V 7] n S" = {(5, 7}. □ 

3.2. AfRne differential contractions and multivariate compositions. For 

the proof of Theorem 11.81 we need to establish first Theorem 13.111 below, which is 
the main purpose of this section. 

The proof of Theorem 13.111 relies on some of the results obtained in [35]. Let us 
introduce the following notation. Given a, 6 G C, 1 < i < j < n and 

F{zi, . . . , Zn) = y] aaZ°' e C[Z1, . . . ,Zn] 



ctSF 



let 



F[ zi, . . . , Zi-i, azi + b-^, Zi+i, . . . , Zj, . . . , Zn ] (5) 



denote the polynomial 



J2 aaz^^ ■ ■ ■ zfH' [azi + b— 



The next lemma follows from [39, Lemma 2.3] by a rotation of the variables. 

Lemma 3.8 (([39J)). If Pq{v) , Pi{v) e C[w] with Pq{v) + xPi{v) 7^ /or 3m(z;) > c 
and '3m{x) > d then 

Po{v)+(x-^\pi{v)^Q 



for 3m(w) > c and '3m{x) > d. 



dv 



Using Lemma \3M and the Grace- Walsh-Szego Coincidence Theorem [51] one can 
argue as in the proof of [33 Proposition 2.2] to show: 

Proposition 3.9 (CW)). Let (ci, . . . , c„) G M" and F e C[zi, . . . , z„] be such that 

F{zi,...,Zn)^0 

if Jm(zfe) > Cfe, 1 < k < n. Then for any 1 < i < j < n one has 
f d \ 

V dzj J 
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whenever '3va{zk) > Ck, 1 < k < n. 

From Proposition 13.91 we immediately get the following. 

Corollary 3.10 (( 39 )). Suppose that 1 < i < j < n and F{zi, . . . , z„) G H„(C). 
Then F{zi, . . . -dj,Zi+i, ...,Zj,..., z„) € 7^„_i(C) U {0}. 

We can now prove the following extension of a famous composition theorem 
of Schur [51] and related results of Malo-Szego [T71 IM] to the multivariate case. 
Further consequences of Theorem 13 . 1 1 1 will be given in Section [4731 and Section [44l 

Theorem 3.11. Assume that all the zeros of f{z) = Yll^QCiiz'' € 7ii(M) are non- 
positive and that F(zi,...,z„) = X]j=o '3i(-^2, . ■ • , ^^n)-?! G Hn(C) and set m = 
min(r, ,s). Then 

ra 

k\akQk{z2, . . . , Zn)z'l e H„(C) U {0}. 

k=0 

Proof. Suppose that / has only nonpositive zeros. Then /(— zqwo) € 7^2(18.) so that 

G{wo,zo, . . . ,z„) := f{~zoWQ)F{zi, . . . ,z„) e 7Y„+2(C). 
By CoroUarv 13.101 we have that 

H{zq, Zi, . . . , Zn) : 

This means that 

H{zi,0, Z2,..., Zn) 

as required. □ 

3.3. Proofs of Theorems 11.81 and II. 9i We can now settle the classification of 
multivariate multiplier sequences stated in Theorem 11.81 

of Theorem \1.8[ For the "only if" direction what remains to be proven is the state- 
ment about the signs. If it were false for some (multivariate) multiplier sequence A 
then since A is a product of univariate multiplier sequences whose entries either all 
have the same sign or alternate in sign there would exist a € N" such that X{a) ^ 
and A(a + ei)X{a + ej) < 0. Let T be the corresponding (diagonal) operator and 
apply it to z"(l — ZiZj) e H„(]R). By Lemma [531 we get 

T(."(l - ...,)) = A(a)z" (l - ^^^^±^11^....) . 

Since 1 + aziZ2 G 7^2 (R) if and only if a < this is a contradiction. 

Now a I— > A(q;) is a multiplier sequence if and only if a i-^- (— l)'"'A(a) is a 
multiplier sequence so we may assume that A (a) > for all a E W\ By applying 
the Ai's one at a time we may further assume that A^ = 1 for 2 < i < n. Hence we 
have to show that if f{zi, . . . , z„) := J^ito Qi{z2, ■ ■ ■ , Zn)z\ € H„(M) and A : N ^ R 
is a nonnegative univariate multiplier sequence then 

M 

. . . , zn)z\ e W„(M) U {0}. 

i=0 



G[ - — ,zo,...,z„ ) G?^„+i(C)U{0}. 



k=0 



klakQk{z2, . . . , zn)z^ e HniC) U {0}, 
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By Theorem 11.71 there are polynomials 



i=0 



with only nonpositive zeros and 



lim pk{z) = V — — 



Z 

1=0 



Furthermore, by Theorem 13. Ill we know that 

M 

J2 >^^MZ2, Zn)z\ e HniR) U {0} 

i=Q 

for each k. Since limfe^oo >^i,k = A(i) we get Y^fLo K'i)Qi{z2, ■ • • , Zn)zl e 7Y„(M) U 
{0} by Lemma [^751 which settles the theorem. □ 

Let us finally prove the characterization of finite order multiplier sequences. 

of Theorem \1.9l The "if" direction is an immediate consequence of Theorem 11.31 
the proof of which is given in Section[4]below. To prove the converse statement note 
first that Theorem 11.81 implies that the symbol of T is the product of the symbols 
of the corresponding univariate operators. Hence it suffices to settle the case n = 1. 
Let F{z,w) — Y^^=o'^kZ^'w^ be the symbol of T. By Theorem 11.71 we know that 
all zeros of 

M 

gm{z) = T[(l + zr] = ak{m)kz''il + 

k=Q 

are real and have the same sign. Note that these zeros are actually nonpositive 
since z = — 1 is a zero of gm{z) for all large m. Now 

M , . 

g„,{z/m) = (l + -)™ E + ^M"' 

V m/ ^ — ' m"^ 

■^^T^ = 1 we have 

M 



lim gm{z/m) = e^y^ Ukz''. 



k=0 

Hence by Lemma 12.31 the polynomial X^feLo'^fc^'^ nonpositive zeros and the 

theorem follows. □ 

4. Algebraic and geometric properties of stability preservers 

4.1. Sufficiency in Theorems [TT2] and [TT3l Since H„(M) C H„(C) it is enough 
to prove only the sufficiency in Theorem 1 1.21 Recall the "affine differential contrac- 
tion" of a polynomial g C[zi, . . . , z„] defined in (O and note that the following 
consequence of Corollarv l3.10l actuallv settles the sufficiency part in Theorem 1 1.21 

Corollary 4.1. Let T e ^n[C] and suppose that Ft{z,—w) G H2n(C). Then 
TeAniC). 
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Proof. If FTiz,-w) e H2„(C) and f{v) G H„(C) then FT{z,-w)f{v) e H3„(C). 
By Corollary 13.101 if we exchange the variables WiS for —d/dvi's the resulting 
polynomial will be in 7i2n(C) U {0}. If we then replace each variable Vi with Zi 
for all 1 < i < n, we get a polynomial in 7i„(C) U {0}. This polynomial is indeed 
T{f). □ 

4.2. Necessity in Theorems [TT2l and [TTSl Let T Ea,/3 aa/s-z"^'^ G X[M]. We 
may write T as a finite sum T = X/7eZ" •^''-^7' where = g^+g.gz ^a''. It 
follows that acts on monomials as T^{z°') — A-y(a)z" for some function A-y : 
N" — > R. The following lemma gives a sufficient condition for A-^ to be a multiplier 
sequence. 

Lemma 4.2. Let T = z'^'T-y G -4„(M) and denote by CH{T) the convex hull of 
the set {7 : T-y ^ 0} . If n £ Z" is a vertex (face of dimension 0) of CH{T) then 
is a multiplier sequence. 

Proof Let v G K!^. If /(z) G H„(R) then f{vz) = /(wiZi, . . . , t;„z„) G 7^„(K) and 

a,/3 

Hence 

:= ^ aa/3w""'^z"9'^ = ^v'^z^'T^ G XW- 

a,/3 7 

Let (z, /i) = a be a supporting hyperplane of the vertex k. Hence, up to replacing 
/X with —fi, if necessary, we have (7 — k, ^) < for all 7 G CH{T) \ {n}. Now let 
V, = i;,(t) = e^**. Then 

7#K 

By letting t 00 we have that zT^ G .4„(M) and the lemma follows. □ 

Let / — E^ a{a)z°' G M[zi, . . . , z„]. Define the support supp(/) of / to be the 
set {a G N" : a{a) ^ 0} and let d — max{|Q;| : a G supp(/)}. We further define 
the leading part of / to be a{a)z'^ , where a is the maximal element of the set 
{7 G supp(/) : I7I — d] with respect to the lexicographic order on Z". Similarly, 
a T — E^ z'l'T^ G .4„[R] let k = maxllaj : ^ 0} and let kq be the maximal 
element of the set {a : \a\ = k,Ta ^ 0} with respect to the lexicographical order. 
Since kq is a vertex of CH(T) we know that A^^ is a multiplier sequence with a 
finite symbol whenever T G We say that T^o is the dominating part of T. 

Note that the dominating part of fg is the product of the dominating parts of / and 
g. Moreover, if X^^{a) 7^ then the dominating part of T{f) is AK.„(a)a(Q;)z"+''°, 
where a(a)z" is the dominating part of / and is the dominating part of T . 

We are now ready to prove that a real stability preserver also preserves proper 
position. Equivalently, Theorem 14.31 below asserts that .4„(M) C An{'C). 

Theorem 4.3. Suppose that T G yl„(R) and that f,g £ 7Y„(M) are such that 
f « g. Then either T(/) < Tig) or T{f) = T{g) = 0. 

Proof. Let T^g be the dominating part of T. We first assume that / = a^z", g = 
Ea baz" G Hum are such that / « g, < deg(/) < deg(ff) and T^o{f)T^oi9) ^ 0. 
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Let the leading parts of / and g be a(a)z" and b{(3)z^ ^ respectively. Let 
fn ^ ^ Cqz" and qh ^ ^ &q-z"- 

|a|=dcg(/) |Q|=dog(g) 

By Hurwitz' theorem Jh and gu are stable. Moreover all coefficients is fu (and in 
gn) have the same sign by [ini Theorem 6.1]. Consider f{vt),g{vt) G 7ii(M) where 
V S M" . Then deg g{vt) = degf{vt) + 1 and the signs of the leading coefficients of 
g{vt) and f{vt) will be the same as the signs of 6(/3) and a{a), respectively. Since 
also f{vt) ^ g{vt) we infer that a{a)b{(3) > 0. 

Now since T^g{f)Ti^g{g) ^ it follows that the leading parts of T{f) and T{g) 
are A„o(a)a(a)z''°"'"" and A^^ (/3)&(/3)z''°+'', respectively. By Theorem ll.Gl fthe mul- 
tivariate Obreschkoff theorem) we know that either T{f) < T{g) or T{g) < T(/). 
As pointed out in the paragraph preceding Theorem 14.31 dominating parts are nec- 
essarily multivariate multiplier sequences and so by Theorem 11.91 we have that 
XKg{a)XKg{f3) > 0. From the above discussion it follows that for v S M" we have 
T{f)ivt) « Tig){vt) with deg(r(/)(i;t)) < deg(T(5)(z;t)), so that T(/) < Tig). 

If deg(/) > deg{g) we may simply repeat the arguments using — / and g. If 
deg(/) = deg(5) we consider / and g + ezif with e > 0. Indeed, deg(/) < deg(5 -I- 
ezif) and f g + ezif by Lemma We then apply the argument of the first 
case, and obtain the desired conclusion as e 0. 

Suppose now that T^^ (/)'7ko {g) — 0. There is nothing to prove if fg = 0. Let 
h^{zi, ...,Zn) = il + ezi)^^ • • • (1 -I- ez„)^" with 6 G N, 1 < i < n, and let /, = hj 
and g^ = h^g. If ^ = (^i, . . . ,^„) is large enough then T^,„{f^)T^„{g^) ^ 0. The 
theorem follows from Lemma [^751 bv letting e ^ 0. □ 

For the proof of necessity in Theorems 11.21 and 11.31 we need to establish first a 
key property for symbols of (real) stability preservers. 

Lemma 4.4. Suppose that F{z,w) G R[zi, . . . . . . ,Wn] is the symbol of an 

operator in y^„(K) and let A G (0,1)". Then F{z,Xw) is also the symbol of an 
operator in An (R) . 

Proof. Suppose that T e yl„(R) has symbol F{zi, . . . , z„, wi, . . . , Wn)- We claim 
that if (5 > then the linear operator SfT defined by 



(5™zi"r(a{'7) 



m=Q 

is an operator SfT : 7i„(R) 7i„(IR)U{0}. If is the linear operator with symbol 
F{zi, . . . , Zn, wi/(l + (5), . . . , Wn) then a simple calculation shows that 

Ts{f)=£fTif{z,il + S),...,Zn)). 

Hence the claim would prove the lemma. 

In order to prove the remaining claim let (5 > and define a linear operator 
RsT : . . . , z„] M[zi, . . . , z„] by 

RsT{f)^T{f) + 6ziT{dif). 

Suppose that / e 7Y„(M) and that T{dif) ^ 0. Since 1 — iw\ G 7i„(C) we know by 
CoroUaryOthat 1+idi G An{C), so f+idif G H„(C), i.e., dif < /. By Theorem 
SJwe know that T{dif) < T(/) and T{dif) < ziT{dif), which by Lemma [U 
gives i?5T(/) G H„(M) U {0}. If T{dif) = then i?5T(/) = T(/) G H„(R) U {0}, 
so RsT G A»(R). 
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An elementary computation shows that when we apply Rs to T m times we get 

™ / \ 

By induction, i?^-^^ : 7^„(M) ^ 7i„(K) U {0} for aU m e N. Now 

m J \ m J \ m J fc! 

It follows that R^/yyiT{f) tends uniformly to £^T{f) on any compact subset of C" 
asm^oo. Thus £^T : 7^„(M) -> 7^„(K) U {0} by LemmaO □ 

From Lemma 14.41 one can easily see that symbols of (real) stability preservers 
actually satisfy the following homotopical property: 

Theorem 4.5. // F{z^ w) € R[zi, . . . , Zn, wi, . . . , m„] is the symbol of an opera- 
tor in ^„(R) then F{^z,Xw) is also the symbol of an operator in ^„(R) for any 
(/i, A) G [0,1]" X [0,1]". Moreover, the corresponding statement holds for symbols 
of operators in ^„(C). 

We now have all the tools to accomplish the proof of the necessity in Theorem 

[Ql 

of Theorem \1.S[ The final step in the proof is to show that F{z, ^iz~^) ^ whenever 
F is the symbol of an operator T E ^„(R), /i € R" and z — (zi, . . . , z„) e C" is 
such that 3m{zi) > for all 1 < i < n. Indeed, since F{z,w) is the symbol of a 
real stability preserver if and only if F{z + a, w) is the symbol of a real stability 
preserver for all a E R" the claim implies that F{z + a,/iz^^) ^ whenever F 
is the symbol of an operator T E An(S.), a € R", G R" and z G C" is such 
that 3m{zi) > for 1 < i < ?i. But it is straightforward to see that any pair 
Z,W E C" such that 3m{Zi) > and 3m(Wi) < can be written as Zi = ai + Zi 
and Wi = HiZ~^, where 3m{zi) > 0, at E R and /i^ > for all 1 < i < n. Thus, the 
theorem follows from this claim. 

Let T = J2a pO-apz^d^^ E An(]^) and let F be its symbol. By multiplying with 
a large monomial we may assume that a^p = if a ^ /3. Let v E R" and denote 
by the operator with symbol F{z,vw). By Lemma 14.41 we have that 

-T{z^)z-^ = J2^o.pv''z'''^{j)p 

= J2 aMvi fz"~^h)pi~^ e u {0} 

for all V E (0, 1)". Fix /i G R" and let v in the above equation be of the form 
with 7 = (71, ... , 7„) G N", where = {"fi^, . . . , 7^"^)- Then v E (0, 1)" for large 
7. Letting 7 tend to infinity and observing that (7)^7^^^ — > 1 we find by Lemma 
Othat 

J2 aaf}l^^z"-f^ = F{z, M^-i) G 7^„(R) U {0}. 

a,0 
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We have to prove that F{z,fj.z ^) is not identically zero. To do this observe that 
F{z, fiz''^) = 51 (^0+^^,(31^'^ ■ 

By Lemma 1121 the dominating part, T^q = Z]/3 "/3+ko,/3 

is an operator as- 
sociated to a multiplier sequence with finite symbol. Hence the nonzero coefficients 
a^+Ko,/3 are all of the same sign by Theorem II .91 This means that the coefficient of 
^Ka F(^z, ^iz~^) is nonzero and proves the theorem. □ 

The proof of necessity in Theorem If .21 now follows easily. 

of Theorem M.SX Let T e ^„(C) and write the symbol of T as w) = Fff(z, w) + 
iFj{z,w), where Ffi{z,w) and Fi{z,w) have real coefficients. Let further Tr and 
T/ be the corresponding operators. Now T : H„(M) H„(C) U {0} so by Lemma 
Owe have that Tr + Zn+iTi : n„{R) 7i:„+i(K) U {0}. Hence by Lemma O 
we know that Tr + (Azi + a)T/ e yl„(K) U {0} for every A S M+ and a € M. 
Suppose that Tr + (Azi + a)Ti = 0. Then T = (i - a - Xzi)Ti, so T/ = since 
i-a-Xzi^Tii (C). We thus have Tr + (Azi + a)Ti e yl„(R) for every A € M+ and 
a e R, which by Theorem Ogives Fr + (Azi + a)Fj e Hn (R) for every A S R+ and 
a e R. By Lemma O and Lemma O this implies that F = Fr + IFj £ 7i„{C), 
as was to be shown. □ 

4.3. The Weyl product and Schur-Malo-Szego type theorems. The results 
of Section 13.21 provide a unifying framework for most of the classical composition 
theorems for univariate hyperbolic polynomials [41 } [5T | I54j. Moreover, they 
lead to natural multivariate extensions of these composition theorems. Let us for 
instance consider two operators S^T G ^n[C] with symbols Fs{z,w) and Ft{z,w), 
respectively. The well-known product formula in the Weyl algebra [6] asserts that 
the symbol of the composite operator ST is given by 

FsTiz,w)^ V ^d^Fs{z,w)d:FT{z,w). (6) 

kSN" 

This suggests the following definition. 

Definition 4.1. Let z = (zi, . . . , Zn), w = {wi, . . . , The Weyl product of two 
polynomials f{z,w),g{z,w) € C[z,i(;] is given by 

{f*g){z,w)^ J2 ^-^d:f{z,w)dZ9{z,w). 

Theorems lf.2l and lf.3l and ^ imply that the Weyl product of polynomials defined 
above preserves (real) stability: 

Theorem 4.6. If f(z,w) and g{z,w) are (real) stable polynomials in the variables 
zi, . . . , z„ and wi, . . . , w„ then their Weyl product (f -k g)[z, w) is also (real) stable. 

Example 4.1 ((Schur-Malo-Szego theorem)). Suppose that S,T G ^i[R] are such 
that Fs{z,w) = /{X^^zw) and Ft(z,w) — g{Xz), where / € 7ii(R) has only 
nonpositive zeros, g G 7^i(R) and A > 0. Then S,T e Ai(R) by Theorem [Ol 
Hence ST G Ai (R) and therefore 

Fst{z, -w) = J2 ^^V^'^(-A-iz«;).g('=)(Az) e K2(R) 

fc>0 
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by ^ . Letting w = and A ^ it follows that 

fc>0 

which is a well-known result of Schur 54] , Malo and Szego [T71 [5T] . 



From Theorem 14.61 one can also recover de Bruijn's composition results |12[ I13j . 
As for composition (or Hadamard-Schur convolution) theorems in the multivariate 
case, we should point out that in [31] Hinkkanen obtained such a result for multi- 
affine polynomials - i.e., multivariate polynomials of degree at most one in each 
variable - that are nonvanishing when all variables lie in the open unit disk D. Note 
that in the case of the open upper half-plane Theorem l4.6l generalizes Hinkkanen's 
composition theorem to arbitrary (not necessarily multi-affine) stable polynomials. 
In fact, by an appropriate conformal transformation one can obtain an analog of 
Theorem l4.6l for multivariate polynomials of arbitrary degrees that are nonvanishing 
when all variables are in D, thus extending Hinkkanen's convolution theorem. 

Finally, by using Theorems 11.21 and 11.31 we can derive yet another property of 
(real) stability preservers: 

Proposition 4.7. Let T G ^n[C] with Ft{z,w) = J2aeN" Qa{z)w"' , where as 
before z = (zi, . . . , z„) and w = (wi, . . . , t(;„). If T <E -4„(C) (respectively, yl„(M) ) 
then Qa{z) e 7^„(C) U {0} (respectively, 7^„(M) U {0}) for all a G N". 

Proof If T e Ai(C), then Ft(zi, . . . , z„, -wi, . . . , -w„) e H2„{C) by TheoremO 
It follows that for any polynomial P{vi, . . . , Vn) G 7i„(C) one has 

P{vi, . . .,v„)Ft{zi,. . . , z„, -wi, . . . , -Wn) e Tizni'C) 

hence 

■ ■ ■ ' " a^)^^^'"'' ■ ■ ■ ' ■ ■ ■ ' """"^ ^ W2„(c) u {0} 

by Corollarv l3.10l Now the polynomial Pa{vi, . . . ,Vn) ■= f"^ • ■ ■ f"" clearly belongs 
to 7i„(C) for any a — (ai, . . . , a„) G N", so that by the above one has 

a\Qa{z) = 

Pa f -TJ^, ■ • • ,~J^)Ft{zi, ...,Zn, -Wl, -Wn) G 7^„(C)U{0}, 

as required. The case when T e -4„(R) is treated similarly. □ 

4.4. Duality, Polya's curve theorem and generalizations. Let us first estab- 
lish the duality property stated in Theorem II. Ill 

of Theorem [rTl\ By Theorem O we have that T e Ai(C) if and only if 

G{z,w) := Ft{z,-w) G Hn{C). 



But Ft* {z, —w) — G{—w,—z) e Tin{C.) so the desired conclusion follows from 
Theorem ll.2l The same arguments combined with Theorem ll.3l prove the analogous 
statement for yl„(R). □ 



In view of Theorem 11.111 one can both recover known results and deduce new 
ones by a simple dualization procedure, as illustrated in the following examples. 
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Example 4.2 ((Hcrmitc-Poulain- Jensen theorem)). Let p{z) — X]fc=o o-kz'' G R[z]\ 
{0}, T — p{d/dz) — J2k=o o.kd'^/dz'' e -4i[M], and let Tp be the linear operator on 
R[z] defined by Tp{f){z) = p{z)f{z). Then T* = Tp so by Theorem [TTI] one has 
T e ^i(M) if and only if Tp e yli(K), which clearly holds if and only if p e Hi(M). 

Example 4.3. The main result of [T] (Theorem 1.4 in op. cit.) shows that any 
operator in Ai (M) that commutes with the "inverted plane differentiation" operator 
= z^D, where D — d/dz, is of the form aD^ for some fc S N and a S M. 
Therefore, by Theorem II. 11 1 we conclude that any operator in ^i(R) that commutes 
with the operator zD^ = is of the form a{zD'^)^ for some fc € N and a e M. 

As we will now explain, both Theorem 11.31 and Theorem II. Ill admit natural geo- 
metric interpretations that lead to further interesting consequences. For simplicity's 
sake, we will only focus on the case n = 1. 

Definition 4.2. Let f{z,w) e R[2;,w] be a nonzero polynomial in two variables of 
(total) degree d and define the real algebraic curve T / (of degree d) by 

Tf^{{z,w)eR^ : fiz,w)^0}. 

We say that /, or equivalently Tf, has the intersection property (1+) if Tf has d 
real intersection points (counted with multiplicities) with any line in of the form 

w — az + (3, where a > 0, /? G M. 

Similarly, we say that / (or Tf) has the intersection property (I-) if has d real 
intersection points (counted with multiplicities) with any line in M'^ of the form 

w = az + P, where a < 0, /3 G M. 

The symbol curve of an operator T G .Aip] with symbol Ft{z,w) G R[z,w] of 
degree d is the real algebraic curve (of degree d) given by 

Tt = {{z,w) G : Ft{z,w) = 0}. 

From Lemma [^T] and Definition 14.21 we get: 

Corollary 4.8. Let J he a nonzero real polynomial in two variables. Then f G 
7^2 (R) */ and only ifTf has the intersection property (I+). 

Therefore, in the univariate case Theorem 11.31 mav be restated as follows. 

Corollary 4.9. Let T e Ai[M.]. Then T e Ai{M.) if and only if its symbol curve 
Tt has the intersection property (I_ ) . 

As depicted in Figure [5] below. Corollary 14.91 essentially allows one to visualize 
whether an operator T e y^i[M] preserves hyperbolicity by checking whether all 
lines in with negative slope has the required number of intersection points with 
Ft. 

In the same spirit, a simple geometric interpretation and proof of Theorem 1 1.1 II 
for n — 1 is as follows: if T G ^i[R] then Ft-'{z,'w) = Ft{w, z) so Ft* is just the 
reflection of Ft in the main diagonal in the zw-plane (i.e., the line w = z). Since 
the intersection property (I_ ) is clearly invariant under this reflection we conclude 
that Ft and Ft* have the aforementioned property simultaneously. 

These geometric reformulations of Theorem ll.3l and Theorem 11.111 provide a uni- 
fying framework for "curve type theorems" that include and considerably strengthen 
Polya's original result [35] and its various known generalizations [TH] ■ 
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Figure 2. Left picture: the symbol curve of degree c? = 6 of 
an operator in yli(R). Right picture: the symbol curve of degree 
d = 7 of an operator in ^i[IR] for which property fails. 



w 




Figure 3. The symbol curve Ft of degree d = 3 of an operator 
T e Ai (K) and its dual curve Ft* 



Example 4.4. In [JH] Polya proved the following result that he considered as the 
most general theorem on the reality of roots of algebraic equations known at the 
time (1916). 

Theorem 4.10 ((Polya's curve theorem)). Let f{x) he a (nonzero) hyperbolic poly- 
nomial of degree n, and let hQ + hix + . . . + bn+mx"'^"^ , where m > 0, be a hyperbolic 
polynomial with bi > for all < i < n. Set 

Gi{x,y) = bofiy) + b,xf'{y) + b^x^f'iy) + . . . + b,J^^\y) S R[x,2/]. 

Then G has the intersection property (I+). 

Proof. By assumption the polynomial q{x) — X]fc=o" ^k^^ has only real and nonpos- 
itive zeros hence q{—xz) g 7^2 (R) and thus the polynomial in variables x, y, z given 
by q{—xz)f{y) belongs to H3(K). From Corollarv l3.10l we then get q{xDy)f{y) e 
7Y2(IK.), where Dy = d/dy, and the result follows by Corollarv l4.8l □ 
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5. Strict stability and strict real stability preservers 

A natural question in the present context is to characterize aU finite order hnear 
differential operators that preserve strict stability and strict real stability, respec- 
tively. These notions are defined as follows: note first that the set of real stable 
univariate polynomials coincides with the set of hyperbolic univariate polynomi- 
als. Denote by 7iJ(M) the set of all strictly hyperbolic univariate polynomials, i.e., 
polynomials in 7ii(R) with only simple zeros. 

Definition 5.1. A polynomial / £ 7i„(K) is called strictly real stable if f{a + vt) g 
Hl{M.) for any a E M" and v G K". One calls a polynomial g G 7i„(C) strictly 
stable if g{zi, . . . , z„) ^ for all n-tuples (zi, . . . , z„) G C" with 3m{zj) > 0. 

Let H,^ (C) (respectively, be the set of all strictly stable (respectively, 

strictly real stable) polynomials in n variables. Clearly, if n > 2 then 7ifj(R) — 
ni{<C) n R[zi, . . . , z„] while (C) n R[z] = M \ {0}. Denote by A' (C) and A' W 
the submonoids of -4„[C] and .4„[R] consisting of all strict stability and strict 
real stability preservers, respectively, i.e., -4^(C) = {T G yln[C] : T{Hf^{C)) C 
H^(C) U {0}} and ^^(R) = {T G X[K] : T(7^^(R)) C WJR) U {0}}. 

In this section we give necessary and sufficient conditions in order for a linear 
operator to belong to either ^^(C) or yl^(R). 

Theorem 5.1. LetT e An[C]. IfT G Af^{C) then Ft {zi, . . . ,z„,-u;i, . . .,-w„) 7^ 
whenever 3m{zj) > and Jm(wfc) > for all 1 < j, k < n. 

Theorem 5.2. LetT e X[R]. IfT G ^^(R) then Ft {zi, . . . ,z„,-u;i, . . . , -Wn) ^ 
whenever 3m{zj) > and 3m{wk) > for all 1 < j, fc < n. 

To prove Theorems 15.11 and 15.21 we need to establish a multivariate extension of 
the following classical result [45l [51] , compare with Theorems 11.51 and 11.61 

Theorem 5.3 ((Strict Obreschkoff theorem)). Let f,g e R[z]. Then 
{af + Pg : a,/3 G M,a2 +/32 > 0} c ^"(M) 

if and only if f + ig e HUC) org + if G HUC). 

Theorem 5.4. Let f, g £ R[zi, . . . , Zn]- Then 

{af + I3g:a,pe M, + /^^ > 0} C H,",(R) 

if and only if f + ig £ "^^(C) or g + if £ H^(C). 

Proof. This is an immediate consequence of Definition 15.11 and Theorem 15.31 □ 

of Theorem[5li Suppose that T G A!;,{C) and let / G HniC). Then 

feizi, . . . , z„) f{zi + ie, . . . , Zn + ie) £ 7Y,^(C) 

for all e > 0, so T{f^) £ H^iC). Letting e ^ it follows from Hurwitz' theorem 
that T(/) G n„{C) U {0} and thus T £ Ai(C). Now Theorem [L2] imphes that 

Ft{zi, . . . , Z„, -Wi, . . . , -Wn) ■■= ^ Qa{zi, 2„)(-w)" G 7Y2n(C), 

where Qa G C[zi, . . . , z„] are not identically zero only for a finite number of multi- 
indices a £ N". Fix (z°, ...,z°)£ C" with 3m(zf ) > 0, 1 < i < n. Then for any 
e > one has 3iTi(z° + ie) > for all 1 < fc < n. Hence 

^ g„(zO + ie,..., zO + ze)(-u.)" G 7^„(C) 



24 



J. BORCEA AND P. BRANDEN 



and by letting e — > we deduce that 

E 



If EcGN" ■ • • ' = then Qa(zO, . . . , z^) = for all a & W and 

consequently T(/)(z2', . . . , z°) = for all polynomials /, which contradicts the 
assumption that T e ^^(C) and Jm(z°) > for all 1 < i < n. □ 

of Theorem \5JA Suppose that T e y^fj(M). If the symbol is not as in the statement 
of the theorem then by arguing as in the proof of the necessity in Theorem 15.11 we 
see that there exists Z° = {z^, . . . , z°) e C" with Jm(z°) > 0, 1 < z < n, and 

F{Z',^w)^0 (7) 

as a polynomial in w — {wi , . . . , w„). Choose a polynomial / e Ti* (M) of sufficiently 
high degree so that T{af + (3dif) ^ whenever o? ^ 0^ > 0. By Theorem 15.41 
we have T{f) + iT{dif) G 7iJ^(C). This is however a contradiction since by ([7]) we 
have T{f){Z°) + iT{dif){Z°) = 0. □ 

The next two theorems give sufficient conditions for operators in the Weyl algebra 
to be strict stability or strict real stability preserving, respectively. 

Theorem 5.5. Let T e -4„[C]. If Ft{zi, . . . , Zn, —wi, . . . , —Wn) e 'H|„(C) then 

Theorem 5.6. Let T e AJK]- // Ft(zi, . . . , z„, -wi, . . . , G H|„(]R) then 

of Theorem 15.51 Let T e An [C] and suppose that 

Ft^zi, . . . ,z„, -Wi, . . . ,-tu„) 

whenever 3m{zi) > and 3m{'Wj) > for all I < i,j < n. If f{vi, . . . , w„) G 'W^(C) 
then 

^t(21, • • ■ -W^l, ■ ■ • , -Wn)fivi, . . . ,t;„) ^ 

provided that 3m(zi) > 0, 3m{wj) > and 3m(tiA;) > for all I < i, j, k < n. By 
Proposition 13 .91 we may replace each variable Wj with Wj — d/dvj for all 1 < j < 
to get 

Ft\Zi,... ,Z„, Wl, . . . , -7\ Wn ] f{vi, . . . ,Vn) ^ 



whenever 3m{zi) > 0, 3m{vi) > 0, 3m{wi) > for 1 < i < n. If we now exchange 
each variable vj by Zj for all 1 < < n, and let for all 1 < i < n, we see 

that 

/ d d \ 

T(/)(zi, . . . ,z„) = FtIzi, . . . , Zn, . . . , j /(zi,...,z„) 7^0 

whenever 3m(zi) > 0, 1 < i < n, hence T(/) e 7i^(C). □ 

of Theorem 15.61 If Ft is as in the statement of the theorem then by Theorem l5.5l we 
have that T E ^^(C). Consider / e 7ifj(M). The case when / is a nonzero constant, 
say f{z) = c e M\{0}, is immediate since T{f){zi, . . . , z„) = cFt{zi, . . . , z„, 0, . . . , 0) ^ 
whenever 3m{zi) > 0, 1 < i < n, hence T{f) e 7i^(M). Suppose that / is not a 
constant polynomial. By re-indexing the variables we may assume that dif ^ 0. 
Then / + idif € 7^^(C), so T{f) + iT{dif) € Hf,(C) U {0}. By Theorem O we 
have that T(/) e H^(R) U {0}, as required. □ 
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To close this section we note that in general the necessary conditions stated 
in Theorems 15.11 and 15.21 are not sufhcient while the sufficient conditions given in 
Theorems I5.5H5.6I are not necessary. This may be seen already in the univari- 
ate case from the following simple examples. The operator T = d/dz is clearly 
strict (real) stability preserving but its symbol Ft{z,w) = w does not satisfy the 
sufficient conditions stated in Theorems 15.51 and 15.61 Consider now the operator 
S — 2z + 1 + [z^ + z)d/dz. One can easily check that Fs{z, —w) ^ whenever 
Jm(z) > and CJm(w) > and also that S preserves strictly real stable (i.e., strictly 
hyperbolic) polynomials. However, 5(1) = 2^; + 1 ^ (C) so S" ^ y^f (C). A char- 
acterization of strict (real) stability preservers would therefore require conditions 
that are intermediate between those of Theorems 15.11 and 15.21 and Theorems 15.51 
and 15.61 



We will now give several examples and applications of the above results. First 
we prove Proposition 1 1 . 1 2l claiming that the polynomial 



with Ai , . . . , An positive semidefinite matrices and B a Hermitian matrix of the 
same order is either real stable or identically zero. 

o/ Provosition ll. 1^ By a standard continuity argument using Hurwitz' theorem it 
suffices to prove the result only in the case when all matrices Ai, A2, . . . , An are 
positive definite. Set z{t) = a + Xt with a G M", A e K" and t e K. Note that 
P := XiAi + . . . + XnAn is positive definite and thus it has a square root. Then 
f{z{t)) = det(P) det{tI+p-^/^HP-^/^), where H := B + aiAi + . . .-|-a„A„. Since 
f{z{t)) is a constant multiple of the characteristic polynomial of the Hermitian 
matrix _ff , it has only real zeros. □ 

6.1. A stable multivariate extension of the Cauchy-Poincare theorem. Let 

A be any n x n complex matrix and define a polynomial C{A, z) = dct(Z — A) G 
C[zi, . . . , z„], where z — (zi, . . . , z„) and Z is the (diagonal) matrix with entries 
Zij — ZiSij. Given I < i, j < n let A^^ be the submatrix of A obtained by deleting 
row i and column j and set Cij{A, z) = det ((Z — A)*-'). For z — {zi, . . . , Zn) and 
I < i < n let z \ z, = {zi, . . . , Zi_i, Zi+i, . . . , Zn), so that Cii{A, z) = C(A", z \ Zi). 

Lemma 6.1. For 1 < j < n one has Tj := 1 + id/dzj £ ^„(C). 

Proof. The symbol Ftj {z, w) of Tj is Ftj (2, — w) = l — iwj and the latter polynomial 
is stable since it is obviously nonvanishing if Jxn{wj) > 0. The assertion now follows 
from Theorem 11.21 □ 

Theorem 6.2. If A is a complex Hermitian nx n matrix then C{A,z) G 7Y„(M) 
and 



6. Multivariate matrix pencils and applications 




C{A",z\zj) < C{A,z) 



for 1 < J < ?i. 
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Proof. Note that since A is Hermitian C{A, z) is real stable by Proposition 11.121 
Now 

C{A^^ ,z\z,) = Cjj {A, z) = z) « z) 

by Lemma [01 and Corollary 12.41 □ 

The above theorem generalizes the classical Cauchy-Poincare theorem stating 
that the eigenvalues of a Hermitian matrix and those of any of its degeneracy one 
principal submatrices interlace. 

An alternative proof of Theorem 16.21 may be obtained by using the following 
consequence of the Christoffel-Darboux identity [77] : 

Lemma 6.3. Let A be any n x n matrix with n>2 and let I < i, j < n. Then 

n 

C{A,y)aj(A,x) ~ C{A,x)C,j{A,y) = Y,iyk ~ Xk)ak{A,x)Ck,(A,y). (8) 

k=l 

Proof. Let X — (xiSij) and Y = (xiSij). The identity 

[X - A)-^ -{Y- A)-^ ^{X - A)-\Y - X){Y - A)-^ 

obtains by multiplying on the left with {X — A) and on the right with (Y — A). 
Taking the ij-th entry on both sides in the above identity and multiplying by 
C{A, x)C{A, y) yields formula ^. □ 

Now let A be a complex Hermitian nx n matrix with n > 2 and let y — z, x = z 
and z = j in dH). Note that Cy(A, z) = Cji{A, z) and since Cii{A, z) = C(A", z\zi) 
we get 

n 

3m{C{A,z)C{A'\7\^)) = Y,^Jm{zk)\C,k{A,z)\^. (9) 

fe=0 

Theorem l6.2l is obviously true for n = 1 and the general case follows by induction 
on n. Indeed, let 3m{zj) > for 1 < j < where n > 2. By the induction 
hypothesis we have C{A^'\ z\zi) £ 7i„_i(C) and then from ^ we deduce that 



^ / C{A,z) \ ^ C{A,z)C{A'\z\zA\ ^ , , „ 
\C{A-\z\z,)) \ \C{A-\z\z,)\^ j- 

Hence C{A, z) and the desired conclusion follows from Lemma [221 

6.2. LzLx conjecture for real stable polynomials in two variables. Here we 
will prove that all real stable polynomials in two variables x, y can be written as 
± det(a;yl + yB + C), where A and B are positive semidefinite (PSD) matrices and 
C is a symmetric matrix. The proof relies on the Lax conjecture that was recently 
settled in [3S] by using in an essential way the results of |5U] . 

Theorem 6.4 (f[30[ [38])). A homogeneous polynomial p g z] is hyperbolic 

of degree d with respect to the vector e = (1,0,0) if and only if there exist two 
symmetric d x d matrices B, C such that 

p(x, y, z) = p{e) det(a;/ + yB + zC). 

We will also need Proposition II . II that we proceed to prove. 
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of Proposition \l.l[ If / G M[zi, . . . , Zn] is of degree d then its homogenization - i.e., 
the unique homogeneous polynomial /h G ^[zi, . ■ . , Zn+i] of degree d such that 
fnizi, . . . , z„, 1) = /(zi, . . . , z„) - is simply 

fnizi, • . • , Zn+l) = Z,t+i/(2l2„+i, • ■ ■ , ZnZ^l{). 

If is hyperbolic with respect to every vector v S R"+^ such that w„+i = and 
Vi > 0, for all 1 < i < n, it follows in particular that for any a — {ai, . . . , an) G K" 
and {vi, . . . ,Vn) G M" the univariate (real) polynomial in t given by 

///((ai, . . . ,q;„, 1) + (wi, . . . ,i;„,0)i) = /(a + 

is not identically zero (since limt_^oc t^'^f{ct + vt) — fnivi, . . . ,Vn,0) ^ 0) and has 
only real zeros. Hence it belongs to 7ii(R). Thus / G 7i„(M) by Lemma \2A\ 
Conversely, suppose that / G Hni^) has degree d and is given by 

/(z) = ^ a„z'', z = {zi, . . . ,Zn). 

Let a = {ai, . . . , an+i) G R"^"'" and v — {vi, . . . , f„+i) G R""'""'" with = and 
Vi > for all 1 < i < n. Since 7^ for some k G N" with |k| = d, Hurwitz' 
theorem yields 

g{z) lim t-''f{tz) ^ V a^z"^ G H„(R). 

kSN", |K|=ci 

Moreover, g is a homogeneous polynomial so by the "same phase property" es- 
tablished in [161 Theorem 6.1] all nonzero a^'s with |k| = d have the same sign. 
Therefore 

fH{v)^g{vi,...,Vn)^ ^ a„uf 1 • • • <" 7^ 

kGN'S \K\=d 

since Vi > for all 1 < i < ^- Now, if q;„+i = then the univariate polynomial 

t ^ fnia + vt) = g{ai + vit, . . . , a„ + tvn) 

has only real zeros by Lemma |2. 11 while if an+i > then again by Lemma |2. li the 
univariate polynomial 

t ^ fnia + vt) = Q;fj+i/(Q;ia,7|i + via'l^^t, . . . , ana~l^ + Vna~l^t) 

has only real zeros. By the last part of Lemma l2.1l the same holds when an+i < 0. 
Hence fn is hyperbolic with respect to all vectors v G R"+^ as above. □ 

Lemma 6.5. Suppose that p G R[a;,?/,2:] is a homogeneous polynomial of degree d 
which is hyperbolic with respect to any (wi,f2,0) G R"^ with vi,V2 G R+. Then 

d 

p{x,y,Q) = ^aix'^-'y\ 

i=0 

where the Oi 's are such that X]f=o ^i^^ polynomial with only nonpositive zeros. 

Proof. By letting z — > it follows from Hurwitz' theorem that p(a;, 0) is hyper- 
bolic with respect to all v G R^, so it is stable by Proposition ll.il By 16, Theorem 
6.1] all the coefficients have the same sign. Since p(x,i/,0) is stable all the zeros 
of the polynomial where p(l, t, 0) = ao -I- ait -I- • • • -I- Odf^ are real and since the 
coefficients have the same sign all the zeros are nonpositive. □ 
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Theorem 6.6. A homogeneous polynomial p G M.[x, y, z] of degree d is hyperbolic 
with respect to all vectors of the form {vi,V2,0) with vi,V2 G ^+ if and only if there 
exist two positive semidefinite dx d matrices A and B and a symmetric dx d matrix 
C such that 

p{x, y,z) — a det(x^ + yB + zC), 
where a G R. Moreover, A and B can he chosen so that A + B = I . 

Proof. Let p be hyperbolic of degree d with respect to all vectors of the form 
{vi,V2,0) with wi, 1)2 G M+ and let a := p{l, 1, 0) G R\{0}. Consider the polynomial 
f{x, y, z) = p(x, X + y,z). Then f{x, y, z) is hyperbolic of degree d with respect to 
all vectors (di,W2,0), where vi,V2 € R+- Moreover, it is hyperbohc with respect 
to the vector e = (1, 0, 0). Hence by Theorem 16.41 there exist two symmetric d x d 
matrices B, C such that 

f{x, y, z) = /(e) det(a;/ + yB + zC). 

Since / is hyperbolic with respect to all vectors of the form {vi , W2 , 0) with vi,V2 G 



K+ we know by Lemma 16.51 that all the eigenvalues of B are nonnegative. Hence 
J5 is a PSD matrix. Let A^ I - B. Then 

p{x, y,z) = a det{xA + y{I — A) + zC), 

and by Lemma 16.51 all zeros of the polynomial 

r(<) := a-^p{l, t, 0) = (1 - tf det (^A + J^I^ ^ R[t] 

are nonpositive. Inverting this we have 

det{A + tI) = (l + trr(^j^^ 

which implies that A has only nonnegative eigenvalues, so that A is a PSD matrix. 

□ 

From Theorem 16.61 and Proposition 11.11 we deduce the following converse to 
Proposition 1 1 . 121 for real stable polynomials in two variables. 

Corollary 6.7. Let f{x,y) G M[a;,?/] be of degree n. Then f is real stable if and 
only if there exist two n x n PSD matrices A, B and a symmetric n x n matrix C 
such that 

f{x, y)=± det{xA + yB + C). 

6.3. Hyperbolicity preservers via determinants and homogenized sym- 
bols. Using Theorem 11.31 with n — 1 and Corollary 16.71 we immediately get the 
following determinantal description of finite order linear preservers of univariate 
real stable (i.e., hyperbolic) polynomials. 

Theorem 6.8. Let T G y^i[K]. Then T G Ai{M) if and only if there exist a G 
d G N, two positive semidefinite d x d matrices A and B and a symmetric d x d 
matrix C such that 

T = adet{zA- wB + C) 

w—d/dz 

From Theorem 16.81 and Proposition ll.il we deduce yet another characterization 
of univariate hyperbolicity preservers involving real homogeneous (Carding) hyper- 
bolic polynomials in 3 variables: 



MULTIVARIATE POLYA-SCHUR CLASSIFICATION PROBLEMS 



29 



Theorem 6.9. Let T G Ai [M] with symbol Ft{z, w) of degree d and let FT{y, z, w) 
he the (unique) homogeneous degree d polynomial such that Ft{1, z,w) = Ft{z,w). 
Then T Cz Ai (M) if and only if the following conditions hold: 

(i) Fr^y, z, w) is hyperbolic with respect to (0, 1, 1), 

(ii) all zeros o/Fr(0,t, 1) lie in (— oo,0]. 
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